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CHAPTER  1 

PROBLEM  DEFINITION 


1.1  INTRODUCTION 

The  objective  of  this  paper  is  to  qualitatively  determine  the  effect  changes  in 
bottom  contour  will  have  on  sound  propagation  using  ray  trace  techniques.  This 
study  will  describe  some  bathymetric  properties  of  the  Nanoose  Bay  range.  It  will 
outline  the  Hamiltonian  Acoustic  Ray  Tracing  Program  for  the  Ocean  (HARPO),  a 
ray  theory  based  model,  and  discuss  the  relationship  between  ray  tracing  and 
transmission  loss.  HARPO  ray  traces  will  be  used  to  portray  the  changes  in  the 
propagation  of  sound  as  changes  in  bottom  contour  are  introduced.  In-situ 
measurements  acquired  from  three  receivers  over  a  relatively  flat  contour  are  also 
examined. 

Staal  and  Deshamais  [1]  used  a  propagation  loss  model,  (PROLOS), 
developed  by  Defense  Research  Establishment  Atlantic  (DREA),  Dartmouth  (Nova 
Scotia)  in  relation  with  in-situ  measurements  in  Nova  Scotia,  over  both  a  smooth  and 
rough  seabed.  PROLOS  is  a  range-dependent  normal-mode  program.  There  are 
many  other  propagation  models  available  that  are  not  range  dependent,  nor  do  they 
address  the  bottom  contour,  some  of  these  arc  listed  in  Appendix  A. 

The  Naval  Undersea  Warfare  Engineering  Station  (NUWES)  uses  the 
Nanoose  Bay  test  range,  and  others  like  it,  to  verify  weapon  radiated  noise  and  test 
weapon  tactics.  It  is  important  to  know  about  the  acoustic  properties  of  Nanoose  Bay 
and  how  the  terrain  will  influence  sound  propagation.  NUWES  currently  models  its 
instrumented  test  sites  assuming  a  flat  bottom  contour.  With  newly  available  ray 
trace  models,  it  is  prudent  to  reevaluate  some  old  and  trusted  ray  tracing  methods 
and  to  incorporate  some  contour  into  the  current  ray  trace  model  for  known  areas  to 
better  predict  sound  propagation.  NUWES  analysts  are  concerned  with 
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measurement  accuracy  and  water  conditions  in  relation  to  transmission  loss  over 
short  ranges  for  evaluating  vehicle  radiated  noise.  During  radiated  noise  data 
acquisition,  hydrophone  saturation  or  inadequate  signal  to  noise  could  contaminate 
radiated  data.  Acoustic  analysis  techniques,  described  by  Keys  [2],  are  also  used  to 
assess  vehicle  and  machinery  performance  by  correlating  vibration  data  to  radiated 
noise  characteristics. 

During  vehicle  tests  on  range  a  directional  source  may  be  aimed  towards  a 
steep  slope  or  hill.  The  change  from  flat  to  sloped  bottom  contour  changes  the 
grazing  angle  which  results  in  an  altered  reflection  coefficient.  Reflection  coefficient 
changes  cause  differences  in  ray  angles,  and  thus  move  focal  point  locations.  Active 
sonar  systems  can  be  attracted  to  these  focal  locations.  Change  in  the  sound 
propagation  path  causes  increased  signal-to-noise,  reverberation  and  focusing  effects. 
High  frequency  sound  is  strongly  affected  by  boundary  scattering  and  microstructure 
variability,  but  low  frequency  sound  is  most  affected  by  sea  floor  properties.  This 
will  be  discussed  in  detail  in  the  following  chapters. 

1.2  DESCRIPTION  OF  NANOOSE  RANGE 

The  Nanoose  range  site,  shown  in  Figure  1  |31,  is  an  instrumented 
underwater  tracking  range  used  to  measure  radiated  noi.se  and  evaluate  active  and 
passive  acoustic  systems.  The  unique  silhouette  of  the  area  is  well  suited  for  use  as 
an  example  to  evaluate  the  influence  of  bottom  contour  on  sound  propagation.  It  is 
situated  on  the  western  half  of  the  Strait  of  Georgia,  offshore  of  Nanoose  Bay  and 
Ballenas  Channel,  and  approximately  18  kilometers  (km)  north  of  Nanaimo,  B.C.. 
Canada.  The  maximum  dimensions  of  the  range  are  25.5  km  (east-west)  by  12.4  km 
(north-south).  It  includes  the  north  end  of  a  deep  trough  within  the  Strait  of  Georgia, 
which  extends  from  the  Ballenas  Islands  in  the  northwestern  portion  of  the  range  in 
an  east  by  southeast  direction  to  central  portions  of  the  Strait  located  between  Sand 
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Head  on  the  Fraser  River  Estuary  and  Valdes  Island.  The  sides  of  the  basin  are 
irregular  with  occasional  outcroppings.  The  lower  slopes  (200  -  4(X)  meters  deep) 
have  more  than  10  meters  of  mud.  Bottom  samples  show  40-70%  clay,  5%  sand 
(diameter  0.625mm-2mm),  and  diatomacious  ooze  and  mud  make  up  the  rest.  The 
average  porosity  is  75%.  Middleton  [4],  and  Herlinveaux  Garrison  [5]  discuss  the 
bathymetry,  geography  and  climatology  of  the  Nanoose  range  in  more  detail.  The 
Nanoose  Bay  test  range  was  selected  for  this  study  because  of  its  steep  sides  and 
irregular  terrain.  NUWES  currently  uses  this  range  daily  for  measuring  radiated 
noise  from  anti-submarine  warfare  weapons.  The  transmission  loss  equations 
currently  used  by  NUWES  do  not  take  this  terrain  into  account,  and  thus  may  be  in 
error. 

The  main  loss  mechanisms  are  volume  absorption,  geometrical  spreading  loss 
(spherical,  cylindrical  etc.),  bottom  reflection  loss  and  surface  and  bottom  scattering 
loss.  Volume  absorption  in  sea  water  is  caused  by  viscous  effects  and  chemical 
relaxation,  which  increase  with  increasing  frequency  [6).  The  bottom  is  acoustically 
absorptive  compared  to  the  surface  (this  is  often  referred  to  as  a  "lossy"  bottom). 
There  are  some  acoustic  problems  where  knowledge  of  sound  propagation  as  it  is 
affected  by  bottom  contour  would  be  useful. 

1.3  RAY  THEORY  AND  TRANSMISvSION  LOSS 

The  Eikonal  equation  (1)  is  the  basis  of  ray  theory.  First,  assume  an 
approximate  solution  to  the  wave  equation  in  terms  of  an  asymptotic  series.  Next 
substitute  the  asymptotic  series  into  the  wave  equation.  This  leads  to 

IAW|2  =  n2  =  Co2/c(R,Zi)2  (1) 

where 

AW  is  the  gradient  measured  in  meters 
n  is  the  index  of  refraction 


CQ  is  the  initial  speed  of  sound 

c(R,Zi)  is  the  speed  of  sound  at  .  ange,  R  and  depth  Zj 
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which  determines  the  ray  paths,  and  energy  transportation.  The  change  in  phase  is 
directly  proportional  to  sound  speed  [7]  and  it  is  not  a  function  of  time  or  frequency. 
Along  the  ray  path,  loss  varies  with  geometric  spreading. 

Limit  the  Eikonal  equation  for  the  case  of  the  ray  path  in  the  vertical 
gradient  where  the  water  is  horizontally  stratified.  If  there  is  an  upward  slope,  the 

shortened  depth  would  interact  with  the  rays  causing  increased  bottom  loss  and 
steeper  ray  angles.  If  the  projector  is  located  at  a  depth,  Zj  where  c(Zj)  =  Cq,  n  is 
the  refractive  index  and  Sq  is  the  direction  of  the  ray  at  this  point,  then 

cos9  /  cosGq  =  c/Cq  =  1/n  (2) 

which  is  Snell's  Law. 

Ray  theory  can  be  used  to  solve  the  Helmholtz  equation.  In  ray  theory, 
waves  travel  along  trajectories  called  rays.  The  wave  speed  is  determined  by  the 
medium's  refractive  index.  Rays  bend  in  response  to  gradients  in  the  refractive  index. 
A  ray  entering  a  layer  of  higher  sound  velocity  is  bent  away  from  me  layer,  while  a 
ray  entering  a  lower  sound  velocity  layer  is  bent  into  the  layer.  Ray  tracing  can  be 
either  graphical  or  numerical  for  solving  the  problem  of  computing  ray  trajectories 
through  a  known  spatial  distribution  of  refractive  indices.  For  this  study,  assume 
rays  striking  the  bottom  are  reflected  specularly.  The  velocity  profile  used  in  the 
HARPO  model  was  a  simplified  version  of  the  in-situ  data  measured  on  12  June 
1990.  Seasonal  variations  of  the  sound  velocity  profile  on  the  Nanoo.se  Bay  range 
measured  over  a  year  and  the  12  June  1990  velocity  profile  are  shown  in  Figure  2. 
Variational  changes  in  water  column  parameters  are  discussed  in  detail  by  Helton  [8] 
Seasonal  changes  in  wind,  mountain  river  runoff  and  convective  mixing  caused 
abrupt  changes  in  the  velocity  profiles. 
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Several  assumptions  are  used  to  calculate  the  sound  intensity  along  the 
ray  path.  The  velocity  changes  linearly  with  depth  (ie.,  the  water  is  horizontally 
stratified),  and  sound  velocity  at  the  operational  depth  does  not  vary  more  than  5 
percent  from  the  surface  velocity.  Therefore,  those  rays  which  leave  the  projector  at 
a  moderate  angle  don't  get  steep,  and  the  sine  of  the  angle  can  be  approximated  as  the 
angle  itself. 


SOllfJD  SPEED  (n>/i| 


FIGURE  2.  Typical  Seasonal  variations  of  the  velocity  profiles  at  the 
Nanoose  Bay  Range  [9,  Figure  A.l. ;  10]. 


The  transmission  loss,  TL,  is  given  by 
TL  =  10  log  (SL/ALfop) 


O) 


where 
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SL  =  source  strength  level,  (intensity) 

ALfop  =  measured  level  at  the  face-of-phone,  (intensity). 

The  transmission  anomaly.  A,  is  how  the  transmission  loss  differs  from  spherical 
spreading.  It  is  generally  defined  as 

A  =  10  log  (SL/ALfop*R2),  (4) 

where 

SL/ALfop,  describes  the  spreading  and  is  the  ratio  of  the  source  level  to  the 

receive  level  at  the  face  of  hydrophone 
and 

R  is  the  slant  range. 

The  transmission  loss  along  the  ray  path  is  directly  proportional  to  the  square 
of  the  path  length  if  only  the  spherical  sp»'eading  loss  is  taken  into  account.  The 
intensity  along  a  ray  path  changes  when  the  sound  wave  reflects  from  the  bottom  or 
the  surface  boundary.  The  intensity  also  depends  on  the  sound  wavelength,  the 
amount  of  energy  focused,  and  how  nearly  perfect  the  focusing  is.  This  is  similar  to 
the  focusing  of  light  by  a  lens.  Sound  has  some  travel  time  from  the  transducer  to  the 
location  where  the  rays  cross  at  the  focal  point  (or  caustic  region).  If  the  travel  times 
are  very  different,  pulsed  tones  would  not  interfere  with  each  other  and  there  would 
be  no  sound  focused.  The  intensity  at  the  focal  location  in  a  convergence  zone  is 
larger  than  the  signal  for  a  spherically  spreading  wave  front  at  the  same  range.  For 
this  study,  assume  specular  reflection  (the  ray  hits  the  bottom  at  angle  0  and  reflects 
off  at  angle  -0)  and  no  reflection  loss.  Since  the  angle  of  incidence  is  equal  to  the 
angle  of  reflection,  the  reflected  wave  front  has  the  same  curvature  and  appears  to 
come  from  an  image  source  below  the  interface,  similar  to  the  Lloyd's  mirror  effect. 

There  is  usually  more  than  one  propagation  path  between  a  source  and 
receiver.  Often  one  path  will  dominate,  and  its  TL  will  be  minimum  compared  to 
other  paths. 
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Clay  and  Medwin  [11]  explain  that  when  sound  rays  converge  or  diverge,  the 
change  of  separation  of  closely  spaced  rays  can  be  used  to  estimate  the  sound 
transmission  loss.  A  caustic  or  focal  point  is  the  envelope  of  the  family  of  rays  when 
they  converge.  The  geometry  of  two  rays  separated  by  a  small  amount  but  not 
converging  is  shown  in  Figure  3  [12].  The  transmission  loss  between  the  source 
and  receiver  which  are  separated  by  a  range  R,  can  be  deduced  from  the  ray  trace  by 
determining  the  vertical  separation  of  closely  spaced  rays.  Ah. 

The  power  radiated  at  the  source  within  the  angle  A0  of  two  closely  spaced 
rays  is,  by  definition,  the  same  amount  of  power  which  radiates  at  the  receiver  within 
the  same  rays.  The  transmission  loss  from  the  source  to  the  receiver  is  therefore 

101og(Ii/l2)  =  I01og(AA2/AAi)  (5) 

where 

11  is  the  intensity  at  the  source 

12  is  the  intensity  at  the  receiver 

AA I  is  the  area  in  a  one  meter  sphere 

AA2  is  the  area  included  by  the  ray  pair  at  the  receiver. 

By  substituting  and  applying  Snell's  Law  equation  (5)  becomes 

TL  =  10  log  (Rh  Ah/A0  cos02/cosei)  =  10  log  (Rj,  Ah/A0  C1/C2)  (6) 

where 

Rh  is  the  horizontal  range  from  the  source  to  the  receiver  (m) 

Ah  is  the  vertical  separation  (m) 

A0  is  the  initial  angle  between  the  two  rays  at  the  source  (rad) 

®1  is  the  angle  from  the  horizontal  at  the  source  (rad) 
is  the  angle  from  the  horizontal  at  the  receiver  (rad) 

The  bottom  can  be  a  relatively  lossy  medium,  but  if  the  receiver  and 
tran.sducer  are  close  to  the  bottom,  and  if  there  is  a  fairly  strong  shallow  water  layer 


or  sound  channel  above  them,  the  bottom  can  be  expected  to  play  a  large  role  in  the 
direction  of  the  ray  and  amount  of  total  transmission  loss. 
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FIGURE  3.  A  diagram  of  the  closely  spaced  rays  used  for 
transmission  loss  calculations  [12,  p.l27]. 

Source  level,  SL,  a  function  of  frequency,  is  equal  to  the  absolute  level,  AL, 

at  the  face-of-phonc,  fop,  minus  the  tran'sniissmn  loss,  TL,  which  is  a  function  of 

range  and  frequency.  The  transmission  !(><;<;  is  3  combination  of  the  spreading  loss 

and  frequency  dependent  absorptive  losses,  a.  inherent  in  the  medium  over  some 

range,  R. 

SL  =  ALfop  +  TL  (7) 

TL  =  SprL  +  aR  (8) 

Transmission  loss  is  used  to  calculate  the  frec-field  sound  pressure  level,  SPL. 
The  Bolt,  Beranek  and  Newman  study  of  N68  [b|,  suggested  the  transmission  loss 
for  Nanoose.TLn,  would  be 


TLn  =  19  log  R  +  1  dB  ,  :n<R<:(X)0  yards 


(9). 
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The  equation  currently  used  for  the  Nanoose  range  [13],  takes  into  account 
the  depths  of  the  source  and  receiver.  The  equation  is 

SprL(R.Z)=  10Iog(l/R2  +  l/(R2+4ZsZr))  (10) 

where 

SprL  is  the  spreading  loss  at  slant  range,  R,  and  depth  Z 
Zg  is  the  depth  of  the  source 
Zj-  is  the  depth  of  the  receiver. 

Derivation  of  equation  (10)  is  given  in  Appendix  B.  This  takes  into  account  the 
addition  of  the  signal  along  the  surface  path  only.  The  pressure  from  the  surface  path 
component  will  add  to  the  direct  path  component  so  the  transmission  loss  will  be  less 
than  that  caused  by  purely  spherical  spreading.  Bottom  contour  is  not  accounted  for 
in  either  of  these  formulas. 

The  NUWES  Underwater  Noi.se  Analysis  Facility  (UNAFAC)  uses  the 
following  frequency  dependent  absorption  formula  based  on  the  discussions  of  Clay 

and  Medwin  [11]  and  Urick  [12].  The  absorption  formula  is 

ai  =(43.6*f2/(f2  +  6400)  +  2.55*1 0-^*f2)/  1000  (11) 

where 

a  I  is  the  absorption  coefficient,  (dB/yd) 
and 

f  is  the  acoustic  frequency,  (kHz) 

The  HARPO  subroutine  SLLOSS  computes  absorption  using  the  Skretting-Leroy 
absorption  formula  [14]  as 

a2  =  a*  (o2/ti)i2  +  b  *  a)2/(o>2^  (d^)  (12) 

where 

a  is  0.006  dBAm 

b  is  0.2635  dB/km 

(Oj  is  1000  Hz  and  0>2  is  1700  Hz 
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02  is  the  absorption  coefficient,  (dB/km) 
f  is  the  acoustic  frequency,  (kHz). 

The  difference  in  the  absorption  values  between  equations  (11)  and  (12)  is 
shown  in  Table  One.  The  loss  values  are  close  but  do  not  agree  and  will  not  be 

compared  in  the  qualitative  analysis.  The  HARPHD  user  may  modify  equation  (12) 
by  changing  the  values  of  Q)i  and  0)2  to  find  a  closer  fit,  if  quantitative  absorption 

values  are  required.  These  absorption  loss  equations  are  valid  for  wavelengths 
much  less  than  the  depth  of  the  water  column.  Other  absorptive  losses  include 
attenuation,  boundary  zone,  and  refractive  scattering  losses.  Attenuation  is 
important  at  higher  frequencies.  Surface  boundary  zone  losses  are  affected  by  wind, 
and  refractive  scattering  loss  can  be  caused  by  fine-scale  changes  in  the  sound  speed 
profile.  Only  spreading  loss  will  be  addressed  in  this  paper. 


TABLE  ONE 


Frequency  (kHz) 

«  NUWES  (dB/m) 

«  SLLOSS  (dB/m) 

Difference  (dB/m) 

.1 

7.7  *  10-8 

9.688  *10-7 

8.9*10-7 

.2 

3.06  *  10-7 

3.838  *  10-6 

3.52  *  10-6 

1 

7.7  *  10-^ 

7.375  *  10- 5 

6.61  *  10-5 

2 

3.09  *  10-5 

1.77  *  10-4 

1.46*  10-4 

5 

1.9  *  10-4 

3.862  *  10-4 

1.96*  10-4 

10 

7.6  *  10-4 

8.563  *  10-4 

9.6  *  10-5 

L3.1  Coherence 

Coherency  is  the  degree  to  which  the  signal  remains  unchanged  or 
undistorted  by  the  process  of  transmission  along  some  path,  including  reflection  at 
the  boundary.  Whenever  an  omni-directional  sound  source  is  present,  both  surface 


12 


and  bottom  reflected  waves  may  arrive  at  a  given  point  to  combine  with  the  direct 
wave.  Depending  on  their  relative  phase,  these  waves  may  either  reinforce  to 
produce  greater  pressure  than  the  direct  wave  alone,  or  partially  cancel  to  produce 
less  pressure.  Incoherence  most  often  occurs  with  high  frequencies,  longer  ranges, 
inhomogeneous  water,  and  rough  boundaries.  If  the  surface  is  relatively  smooth 

compared  to  wavelength,  the  reflected  wave  acts  as  if  it  were  emitted  by  an  image 
located  a  distance  2^  above  or  below  the  surface  and  1 80®  out  of  phase  with  the 
source,  where  is  the  depth  of  the  receiver.  An  incoherent  wave  has  a  random 
phase  that  is  uniformly  distributed  over  a  length  period  of  2ji;  its  mean  power 
densities  may  be  added.  A  coherent  wave  has  a  constant  phase  and  a  total  power 
density,  obtained  by  summing  the  individual  fields  vectorially  to  determine  the  power 
density  from  the  resulting  total  field.  For  the  Nanoose  example,  the  pressure 
amplitudes  should  be  added  incoherently. 

Scattering  strength  increases  with  increasing  frequency  and  depends  on 
bottom  and  subbottom  composition  and  roughness.  In  Physics  of  Sound  in  ilis  ^£3 
[15],  scattering  strength  is  described  in  detail.  As  grazing  angles  tend  to  0®, 
scattering  strength  decreases  rapidly  to  become  negligible.  If  the  bottom  has 
irregular  features  such  as  rocks  and  boulders,  one  can  expect  a  large  amount  of 
backscatter.  Clay  and  Medwin  [11]  suggest  that  for  a  1  kHz  sonar  signal,  an 
irregularity  as  small  as  1  meter  is  large  enough  to  strongly  scatter  the  sound. 

Although  the  surface  boundary  is  considered  to  be  where  the  primary  effects  of 
multipath  interference  originate,  acoustics  involving  the  bottom  boundary  is  an 
important  source  of  specular  reflection,  back  .scattering,  absorption  into,  and 
transmission  from  bottom  sediments.  Exact  solutions  require  completely  specified 
geometry  and  reflection  coefficients.  The  refraction  and  transmission  of  sound  at  the 
water/sediment  boundary  can  be  expressed  in  terms  of  the  ratio  of  the  acoustic 
impedances,  Pici/p2C2,  and  angle  of  incidence  with  the  bottom.  After  incidence  on 
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a  rough  boundary,  the  reflected  signal  will  be  reduced  in  amplitude  and  changed  in 
phase  by  the  reflection  coefficient  of  that  boundary.  If  the  value  of  the  speed  of  sound 
for  the  subbottom  appears  low,  it  may  be  due  to  a  sloping  profile  giving  a  low 
apparent  velocity. 

The  ocean  bottom  has  both  topography  and  reflectivity  (material  property 
such  as  porosity).  Newhall  et.  al.,  maintain  that  reflectivity,  which  is  a  function  of 
the  local  geology,  angle  of  ray  incidence,  and  acoustic  frequency  can  be  described 
by  the  plane  wave  reflection  coefficient  [16].  Based  on  the  grazing  angle  and 
knowing  the  composition  of  the  bottom,  the  HARPO  user  would  look  up  the 
reflection  coefficient  in  a  table.  Bottom  loss  is  usually  measured  by  comparing 
direct  path  to  bottom  bounce  path  after  subtracting  spreading  loss  for  path  length 
differences.  Bottom  loss  can  be  expressed  as  20  log  III ,  where  £  is  the  reflection 

coefficient,  and  is  small  for  small  angles.  The  UNAFAC  equations  and  HARPO  do 
not  incorporate  the  reflection  coefficient  at  this  time.  More  discussion  is  beyond  the 
scope  of  this  paper  and  readers  are  directed  to  Morris  et.  al.,  who  explain  that  some 
models  take  into  account  speed  of  compressional  and  shear  waves  in  the  sediment, 
absorption  coefficients,  density,  and  the  thickness  of  the  subbotiom  layers  (17). 

1.4  HARPO 
1.4.1  OVERVIEW 

In  the  past,  various  forms  of  the  parabolic  equation  were  considered  the  only 

i 

way  to  handle  change  in  bottom  contour.  Some  models  have  become  available  that 
u.se  ray  tracing  and  take  this  into  account.  This  .section  is  a  brief  overview  of  the 
theory  behind  the  HARPO  model. 

HARPO  stands  for  the  Hamiltonian  Acoustic  Ray-tracing  Program  for  the 
Ocean.  It  is  a  general  purpose  FORTRAN  ray  tracing  program  which  traces  three- 
dimensional  paths  of  acoustic  rays  through  model  oceans.  It  can  be  used  to  take  into 
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account  the  3  dimensionality  of  oceans  and  to  compute  reflections  from  complicated 
bottom  models.  HARPO  calculates  ray  paths  by  numerically  integrating  Hamilton’s 
equations  in  four  dimensions,  three  spatial  and  one  temporal,  for  earth-centered 
spherical- polar  coordinates  p,9,  <}»  and  time  t,  as  well  as  its  first  spatial  and  temporal 
derivatives.  Jones  er.  a/.,  have  shown  [14j  the  solutions  to  the  wave  equation  are 
related  to  paths  that  satisfy  Fermat's  principle,  and  that  Hamilton's  equations  can  be 
integrated  to  construct  such  paths.  In  Cartesian  coordinates,  Hamilton’s  equations 
take  the  form 

dxj/di  =  5H/5ki  ;  dki/di  =  -  5H/5xi  ,  i  =  1  to  3  (13) 

where 

T  is  sometimes  time,but  depends  on  how  the  Hamiltonian,  H,  is  defined 
kj  are  wave  number  components 
Xj  are  coordinates  of  a  point  on  the  ray  path 

To  solve  equation  (13)  for  the  ray  path,  the  user  chooses  initial  values  for  the 
quantities  of  Xj  and  kj  integrating  results  in  six  differential  equations.  For  acoustic 

waves  in  the  ocean,  the  Hamiltonian,  constant  along  a  ray  path,  is  defined  as 

H(xi,ki)  =  (0)  -  k  *  V(xi))2  -  C2(xi)k2  =  0  (14) 

where 

V(xj)  is  the  ocean  current,  may  be  time  dependent 
C(xj)  is  the  sound  speed  field,  may  be  time  dependent 

and 

0)  is  the  angular  wave  frequency. 

Equation  (14)  shows  that  the  effect  of  a  three  dimensional  vector-current 
field  i?  included  in  the  definition  of  the  Hamiltonian  Hamilton's  equations  are 
equivalent  to  the  Eikonal  equation  for  determining  the  ray  path. 

Hamiltonian  ray  tracing  is  an  alternative  to  .segmented  methods  and  requires 
the  ocean  to  be  modeled  as  a  continuous  three-dimensional  function  with  continuous 
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gradients.  Each  ray  path  is  computed  by  numerically  integrating  Hamilton’s 
equations  with  a  different  set  of  initial  conditions.  Because  it  uses  continuous  models, 
this  method  avoids  false  caustics  and  discontinuous  ray  paths.  These  properties  are 
encountered  in  conventional  ray  tracing  methods,  which  use  layers  where  each 
acoustic  ray  path  segment  is  computed  alone.  The  reader  is  directed  to  Jones,  Riley 
and  Georges  [14],  who  have  documented  the  software.  The  user  defines  an  ocean 
model  with  a  unique  velocity  profile  and  bottom  depth. 

A  version  of  HARPO  has  recently  been  developed  by  Georges  et.  al.  [18]  to 
run  on  a  personal  computer  (PC)  using  a  FORTRAN  compiler  as  an  improvement  to 
the  Digital  Equipment  Corporation,  VAX  version.  The  HARPO  model  will  give  the 
user  freedom  to  manipulate  the  various  parameters  which  affect  the  sound 
propagation,  including  bottom  contour,  current  and  wind  influence.  It  can  display 
ray  traces  both  graphically  and  in  tabular  form. 

HARPO  can  describe  the  geometry  of  the  ray  path,  calculate  a  pulse  travel 
time  and  phase  time,  and  compute  the  Doppler  shift,  rate  of  change  of  phase  (if  the 
user  has  chosen  the  ocean  to  vary  in  time),  absorption  and  geometric  path  length.  It 
allows  continuous  three-dimensional  models  of  the  refractive-index  field  and  two- 
dimensional  models  of  reflecting  surfaces.  The  user  may  trade  computing  speed  for 
accuracy  by  specifying  the  maximum  allowable  integration  error  per  step.  I  have 
only  manipulated  the  bottom  topography  by  introducing  a  Lorentzian  shaped  ridge  to 
the  flat  bottom  model. 

1.4.2  LIMITATIONS  AND  ASSUMPTIONS 

Criteria  for  using  ray  tracing  include  making  sure  there  is  a  fractional  change 
in  gradient,  g,  over  a  wavelength,  X,  which  is  small  compared  to  sound  speed,  c,  (ie.. 
gX  «c)  where  g  =  dc/dz  =  (ci-C2)/(Z|  -Z2),  therefore  ray  models  are  more  accurate 

at  higher  frequencies  (e.g.,  above  1  kHz).  A  second  criterion  states  that  the  radius  of 
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curvature  is  much  greater  than  the  wavelength.  This  means  that  only  plane  waves  in 
the  far  field  are  acceptable.  The  velocity  is  assumed  to  change  linearly  with  depth. 
Because  it  is  based  on  ray  theory,  HARPO  has  the  same  limitations.  The  ocean 
model  the  user  inserts  must  be  deterministic,  not  random.  It  accepts  only  bottom 
models  with  continuous  surface,  slope  and  curvature.  For  example,  wedge  shaped 
surfaces  are  not  allowed.  It  cannot  handle  refraction  at  discontinuities  of  refractive 
indices  or  its  gradients. 

HARPO  does  not  check  to  be  sure  the  input  from  the  user  is  reasonable,  nor 
does  it  apply  corrections  for  diffraction  or  partial  reflections.  It  does  not  check  to  see 
if  ocean  models  satisfy  physical  conservation  laws  and  boundary  conditions.  It  does 
not  check  any  models  to  see  if  they  satisfy  continuity,  or  physical  boundary 
conditions,  nor  does  it  automatically  compute  amplitude  or  transmission  loss. 
Therefore,  the  total  amplitude  at  a  receiver  must  be  computed  by  the  user  by 
combining  the  absorption,  reflection  losses  and  focusing  listed  for  each  ray.  The 
user  must  also  decide  whether  to  add  the  contributions  of  multipath  rays  coherently 
or  incoherently. 

HARPO  assumes  conservation  of  energy  within  the  bundle  of  rays,  so  the 
wave  intensity  is  inversely  proportional  to  the  cross  sectional  area  called  a  flux  tube. 
When  that  area  becomes  zero,  ray  theory  predicts  infinite  energy  density  (a  caustic, 
where  two  or  more  rays  cross). 
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CHAPTER  2 
ANALYSIS 


2.1  HARPO  SETUP 

While  modeling  the  ray  paths,  the  frequency  and  angle  of  incidence  were 
held  constant.  The  difficulty  with  the  topography  was  interpolating  the  data  from 
bathymetry  charts  into  a  form  usable  to  the  code  without  taking  up  too  much  storage 
space  and  thus  increasing  the  program  run  time.  The  rays  of  the  most  interest  are  the 
eigenrays  connecting  the  source  to  the  receiver.  That  is  the  direct  path,  the  one 
surface  bounce,  the  one  bottom  bounce  and  the  one  bottom/surface  bounce. 

"Shooting  techniques"  must  be  used  to  find  these  rays.  The  user  sends  out  a  fan  of 
rays  at  different  launch  angles  which  bracket  the  receiver  range,  depth  and  azimuth. 
These  ray  launch  angles  are  then  linearly  interpolated,  new  rays  are  launched  based 
on  this  first  guess.  This  procedure  is  iterated  until  a  given  ray  has  come  within  a 
specified  distance  of  the  receiver.  The  user  gets  an  accurate  first  guess  for  the 
launch  angles  to  other  ranges  (elements)  by  finding  the  eigenray  for  the  farthest  range 
(element)  radially.  I  was  not  able  to  establish  the  appropriate  eigenrays  to  determine 
the  actual  transmission  loss  levels,  this  results  in  a  strictly  qualitative  analysis. 

The  GLORENZ  subroutine  supenmposes  a  Lorentzian  shaped  ridge  on  a 
spherical  surface  of  any  height  and  width  to  represent  the  steep  slope  and  hills  typical 
on  the  Nanoose  Bay  range.  The  west  end  of  the  Nanoose  Bay  range  was  the 
incentive  for  the  steep  slope  representation  The  bathymetric  profile  of  the  Nanoose 
Bay  range  along  the  centerline  is  shown  in  Figure  4.  The  hilly  bottom  contour  profile 
across  Array  7  is  shown  in  Figure  5.  These  are  models  of  the  most  commonly 
occurring  bottom  type,  beside  the  flat  bottom,  on  the  Nanoose  range  shown  in  Figure 
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FIGURE  4,  The  Bathymetric  profile  of  the  Nanoose  Bay  range 

along  the  centerline. 


FIGURE  5.  The  bathymetric  profile  across  array  7 
on  the  Nanoose  Bay  range. 


Georges  points  out  "the  bathymetry  sensitivity  is  critical  toraypath  stability 
issues,  which  affect  signal  processing,  and  topography"  (18,  p.  203].  There  arc 
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many  rays  both  upward  and  downward  from  the  horizontal,  for  clarity  only  -2  degree 
through  -6  degree  launch  angles  have  been  shown  in  Figures  6  through  8  and  10. 

2.2  HILL  CONTOUR 

The  Lorentzian  shaped  hill  topographies  were  put  into  the  HARPO  model  at 
both  100m  and  200m  high  as  shown  in  the  ray  trace  Figures  6  and  7.  The  larger  hill 
is  similar  to  the  contour  found  on  the  Nanoose  range  in  array  7  ,  as  is  shown  in 
Figure  5.  For  modeling  purposes,  the  area  in  front  and  behind  the  hill  is  flat. 
Compare  the  ray  trace  of  the  flat  bottom  in  Figure  8  with  that  of  the  larger  hill. 

Figure  7.  The  rays  become  steeper  as  they  hit  the  hill.  A  shadow  zone  develops 
beyond  and  behind  the  hill.  The  angle  dependent  bottom  reflection  coefficient 
increases  when  the  grazing  angle  increases.  The  steeper  rays  increase  the  reflection 
coefficient  and  increase  the  bottom  loss  experienced  by  the  sound.  The  ray  paths  are 
warped,  the  path  lengths  are  different,  and  spreading  loss  is  increased.  The  channel 
trapped  ray  is  destroyed. 


FIGURE  6.  A  HARPO  ray  trace  with  small  hill  topography. 
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FIGURE  7.  A  HARPO  ray  trace  with  large  hill  topography. 


FIGURE  8.  A  HARPO  ray  trace  with  flat  topography. 
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The  smaller  hill  does  not  effect  the  ray  as  much  as  the  larger  hill  since  it  is 
below  the  sound  channel.  The  channel  trapped  ray  is  undisturbed.  The  ray  bottom 
interaction  is  dependent  on  the  height  of  the  hill  with  respect  to  the  depth  of  the 
transducer  and  the  shape  of  the  velocity  profile.  Rays  leaving  the  source  at  steep 
angles  interact  with  the  hill  and  suffer  reflection  loss.  This  phenomena  has  been 
studied  by  Jensen,  who  wrote,  Ray  angles  steepen  by  twice  the  bottom  slope  per 
bounce,  ie.  around  28°  and  hence  can  undergo  a  maximum  of  three  upslope 
reflections  before  being  redirected  back  toward  the  source”[19,  p.  193].  The 
example  in  Figure  9  shows  the  near-grazing  rays  get  over  the  top  and  get  trapped  in 
the  upper  .sound  layer.  In  front  of  the  hill,  the  propagation  is  dominated  by  steep 
bottom  bounce  paths.  Pan  of  the  sound  bounces  back  toward  the  source  to  create  a 
large  caustic  or  focal  region.  Volume  reverberation  increases,  causing  a  low  signal- 
to-noise  ratio  due  to  the  rays  bouncing  off  the  slope.  The  volume  of  water  in  front 
of  it  is  ensonified. 


SOUND  SPEED  RANGE  (  km  ) 

(m/a) 


FIGURE  9.  A  ray  trace  of  the  sound  propagation  across  a  seamount 

[18,  p.  193]. 
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Multiphone  averaging  used  to  reduce  the  multipath  effect  will  not  be  enough 
to  reconcile  the  signal,  in  the  case  of  the  larger  hill.  Based  on  the  fact  that  the 
bottom  is  very  lossy,  I  would  expect  the  sound  in  the  flat  contour  to  die  out  more 
quickly.  The  rays  in  the  hill  and  slope  examples  tended  to  get  trapped  in  the  sound 
channel.  In  the  flat  example,  the  rays  just  bend  and  don't  hit  the  bottom,  but  in  the 
hill  example,  the  rays  are  distorted.  The  minus  two  degree  ray,  for  example,  would 
not  hit  the  bottom  with  the  slope,  whereas  on  the  hill  it  does.  Then  the  sound  will 
experience  more  loss  due  to  the  lossy  bottom. 

2.3  SLOPE  CONTOUR 

The  slope  topography  used  in  the  HARPO  model  is  similar  to  the  western 
slope  on  the  centerline  of  the  Nanoose  Bay  range.  The  bottom  contour  of  the  range 
along  the  centerline  is  shown  in  Figure  4.  The  centerline  of  the  range  is  where  most 
of  the  events  tested  on  range  take  place.  Compare  the  flat  contour  ray  trace  in 
Figure  8  to  the  slope  topography  ray  trace.  Figure  10.  The  sound  is  focused  back 
toward  the  source  and  the  volume  of  water  in  front  of  the  slope  is  ensonified.  This  is 
similar  to  the  effect  of  the  hill  but  more  pronounced.  The  rays  angles  are  steepened, 
and  the  sound  has  become  focused. 

Three  Ambient  Noise  Source  (ANS)  hydrophones  (ANS  1,  ANS  2,  and 
ANS  3)  are  permanently  mounted  approximately  12  feet  from  the  floor  of  the 
Nanoose  range  as  shown  in  Figures  4  and  8.  Local  roughness  will  cause  scattering 
because  of  their  close  proximity  to  the  ocean  floor  .  ANS  1  is  nearest  the  base  ot 
the  steep  slope  and  may  be  influenced  by  the  increased  reverberation,  as  the  volume 
of  water  in  front  of  the  slope  is  ensonified.  This  may  cause  a  low  signal  to  noise 
ratio.  As  the  autogain  of  ANS  1  increases  to  compensate  for  the  reverberation,  the 
hydrophone  may  not  detect  signals  or  radiating  noise.  ANS  2  and  ANS  3  are  further 
from  the  steep  slope  and  will  have  a  higher  signal-to-noise  ratio.  In  the  flat  bottom 
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contour  example,  sound  ray  paths  are  considerably  more  influenced  by  the  velocity 
profile  than  the  irregular  bottom  contour  examples. 
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FIGURE  10.  HARPO  ray  trace  with  steep  slope  topography. 
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CHAPTER  3 

TRANSMISSION  LOSS  MEASUREMENTS  USING  ANS 

3.1  MEASUREMENT  SETUP 

The  objective  of  the  test  performed  at  the  Nanoose  Bay  range  on  12  June 
1990  was  to  calibrate  the  three  ANS  hydrophones.  These  hydrophones  were  evenly 
spaced  along  the  bottom  of  the  range  over  a  relatively  flat  contour.  It  was  hoped  that 
these  results  could  also  be  used  as  the  baseline  set  of  measurements  of  spatially 
separated  receivers  over  a  flat  bottom  contour.  This  baseline  can  be  compared  with 
future  in-situ  measurements  with  the  same  receiver  spacing  over  an  irregular  bottom 
contour  to  quantitatively  determine  the  difference  in  transmission  loss  between  the 
flat  and  irregular  bottom  contours.  The  Sonar  Acoustic  Target  Source  (SATS) 
transducer  array  [20]  was  used  to  project  tones  during  the  test.  SATS  is  a  TR- 
292AVQM-6.  SATS  is  a  target  source  and  measuring  system  used  for  the 
certification  of  sonar  systems.  It  is  semi-portable,  surface  deployed  over  the  side 
equipment.  It  is  made  of  several  components,  two  of  which  were  required  to  span  the 
desired  frequencies.  A  random  sample  of  center  frequencies  of  100,  200,  500,  1  kHz, 
2  kHz,  5  kHz,  10  kHz,  and  20  kHz  were  projected  from  SATS.  The  channel  4 
projector,  TR-295/WQM-6,  consisted  of  a  coated  ceramic  sphere  mounted  on  a  base 
plate  through  which  a  cable  and  penetrator  housing  entered.  For  the  lower 
frequencies,  projector  TR-294AVQM-6  was  used.  It  was  composed  of  six  bender  bar 
assemblies  surrounded  by  a  oil-filled  neoprene  boot,  a  compliant  tube  assembly 
surrounded  by  a  rubber  boot,  shields  and  guard  bars,  a  transformer,  and  cover  and 
frame  assemblies.  The  hydrophones  used  to  acquire  the  sound  were  three  ITC- 
H68A  hydrophones.  These  phones  were  on  12  foot  frames  planted  on  the  bottom  of 
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the  Nanoose  Bay,  approximately  6,000  meters  from  each  other  as  shown  in  Figures  1 
and  8. 

During  the  test,  SATS  projected  a  series  of  continuous  wave  signals  from  a 
location  between  ANS  2  and  ANS  3  at  a  depth  of  approximate  250  feet. 

3.2  ANALYSIS  OF  RESULTS 

At  least  10.5  seconds  was  required  for  a  sound  wavefront  to  get  from  SATS 
to  the  slope  and  back  to  ANS  1.  The  analysis  was  averaged  over  a  4  to  23  second 
time  period  during  the  25  to  55  second  broadcast.  One  would  expect  the  level  at 
ANS  1  to  be  different  than  the  other  ANS  hydrophones  if  the  slope  interfered  with 
the  direct  path. 

Table  Two  contains  source  level  estimates  of  the  data  at  the  ANS 
hydrophones.  This  is  the  absolute  level  at  the  face-of-phone  minus  the  Nanoose 
transmission  loss  (equation  10).  There  were  no  low  frequency  results  due  to  the 
surrounding  high  ambient  conditions  caused  by  several  craft  generators.  The  SNR 
was  not  high  enough  to  give  me  confidence  in  the  output  level  at  the  low  frequencies. 
There  were  no  significant  differences  in  the  levels  of  the  three  hydrophones  over  the 
relatively  flat  bottom.  The  level  of  reverberation  was  too  low  to  interfere  with  the 
direct  path.  This  may  have  been  due  to  the  fact  that  the  bottom  and  slope  are 
covered  with  mud  and  ooze,  creating  a  very  lossy  bottom.  No  shadowing  occurred 
because  the  area  between  the  three  phones  was  relatively  flat. 
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TABLE  TW< 

[) 

Freq.  (kHz) 

PAMS 

SLreal  (dB) 

ANS  1 

SL  (dB) 

ANS  2 

SL  (dB) 

ANS  3 

SL  (dB) 

2 

162 

156 

154 

162 

5 

159 

162 

157 

161 

10 

164 

163 

161 

162 

20 

164 

163 

J61 _ 

162 

3.3  RECOMMENDATIONS  FOR  THE  NEXT  EXPERIMENT 

To  measure  the  effects  of  an  irregular  bottom  contour,  an  additional  set  of 
tests  should  include  setting  one  of  the  receivers  on  the  opposite  side  of  the  hill  from 
the  transducer,  where  the  expected  signal  level  will  be  much  lower.  Next,  position 
the  receivers  along  the  base  of  the  steep  slope,  and  place  a  directive  source  beyond  it 
The  directive  source  could  then  be  aimed  in  several  directions  to  measure  the  effect 
of  the  volume  reverberation  and  focusing  off  the  slope.  Another  good  test  scenario 
would  be  to  position  an  omni-directional  source  far  from  the  slope  with  a  directive 
receiver  located  near  the  slope. 
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CHAPTER  4 

CONCLUSIONS  AND  RECOMMENDATIONS 

This  study  focuses  on  ray  trace  deformation  due  to  changing  bottom  contours. 
The  HARPO  program  was  manipulated  to  model  several  bottom  contour  types,  but 
actual  transmission  loss  levels  were  not  obtained.  This  results  in  a  strictly  qualitative 
analysis. 

Changes  in  bottom  contour  have  been  shown  to  affect  sound  propagation 
parameters,  such  as  changes  in  path  length  and  angle  of  incidence.  The  result  is  a 
difference  between  the  predicted  focal  point  (caustic)  locations  for  flat  contour  and 
the  irregular  contour.  The  results  of  in-situ  measurements  are  valuable  as  a 
quantitative  baseline,  and  can  be  compared  to  future  transmission  loss  level 
measurements  acquired  from  receivers  at  the  same  range  separation  over  irregular 
bottom  boundaries. 

Other  variables  can  be  introduced  to  the  HARPO  model,  including  ocean 
currents,  which  can  be  used  to  predict  the  influence  of  the  tidal  currents  prevalent  in 
this  region.  The  reflection  coefficient  for  the  bottom  should  be  added  to  the  program 
capability. 

On  the  Nanoose  Bay  range,  a  directive  beam  in  the  westward  direction  would 
simulate  the  steep  slope  model,  while  the  hill  portrayed  is  a  simulation  of  the  contour 
encountered  in  array  7.  Other  test  ranges  can  be  modeled  in  this  same  manner.  The 
Dabob  range  is  much  smaller  but  has  similar  steep  boundaries.  The  Quinault 
Underwater  Tracking  Range  is  relatively  flat  and  very  shallow. 

Quantitative  transmission  loss  data  will  be  useful  in  further  studies  of  ray 
tracing,  and  would  be  a  valuable  addition  to  the  studies  already  accomplished  on  the 
NUWES  underwater  tracking  ranges. 
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APPENDIX  B 

NANOOSE  TRANSMISSION  LOSS  DERIVATION 


The  Nanoose  range  site  transmission  loss  equation  was  developed  based  on 
geometry  of  a  surface  bounce  adding  with  the  direct  path  [20]  shown  in 
Figure  1 1.  The  surface-bounce  path  range,  Rg  is 

Rs  =  Li-»-L2=  (Rh2-H(Zs+Zr)2)l/2  (Al). 

where  Rj,,  the  horizontal  range  is 

Rh  =  (  r2  +  (Zs  -  Zr)2)  1/2  (A2). 

Then 

Rs  =  (R2  -  (Zs  -  Zr)2  -h  (Zj  -HZr)2)l/2 

=  (R2  -  Zs2  +  2ZsZr  -  Zr2  -h  Zs2  -t-  2ZsZr  -t-Zr2)l/2 
=  (R2  +  4ZsZr)l/2  (A3) 

Assume  only  the  first  order  surface  and  the  direct  signals  come  in  to  the  receiver. 


BOTTOM 


FIGURE  11.  A  Nanoose  transmission  loss  diagram. 
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The  received  signal.  r(t),  then  becomes 

r(t)  =  s(t)/R2  +  s(t)/Rs2  (A4), 

where 

s(t)  is  the  source  signal 
r(t)  is  the  received  signal. 

Assuming  the  surface  is  a  perfect  reflector, 

r(t)  =  s(t)  (  1/r2  +  1/(r2  +  4Z3Zr  ))  (A5). 

The  spreading  loss  becomes 

SprUR.Z)  =  10  log(  1/r2  +  l/(R2+4ZsZr))  (A6) 

where 

SprL  is  the  spreading  loss  at  slant  range,  R,  and  depth  Z 
Zj  is  the  depth  of  the  source 

Zj-  is  the  depth  of  the  receiver. 

Given  the  transmission  loss  is  the  spreading  loss  pluss  the  absorption  loss, 

TL  =  SprL  +  aR 

the  Nanoose  transmission  loss  equation,  TLjsj,  becomes: 

TLn  =  10  log  (1/r2  +  1/(r2  +  4z^Zr  »  +  aR.  (A7) 

The  B,B  &  N  study  of  1968  [9]  suggested  the  transmission  loss  for  Nanoose.TLn, 
would  be 

TLr,  =  19  log  R  +  1  dB  +  aR  ,  20<R<2000  yards  (A8) 

The  difference  between  the  two  equations  begins  to  get  measurable  at  long  distances. 
For  example. 

Given: 

R=  1900  yds 
Zs  =  Zr  =  200  feet 

a  =  8.563  *  10-4dB/yd 

TLjvj  =  67.2  dB  and  TLn  =  64.9  dB,  therefore  ATL  =  2.3 


